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ABSTRACT 


The choice of correctly defined and useful parameters for charac- 
terizing the performance of reference clocks used in digital com- 
munications systems is a subject of efforts in the CCITT and the CCIR, 
the consultative technical bodies of the ITU. The operation of the 
systems depends on timing. Thus the CCITT has chosen in its Recommen- 
dation G811, to specify the clock performance requirement by starting 
a maximum allowable time interval error, i.e. a clocks time departure 
after Initial synchronization. 

On the other hand the commonly used characterization of clock perfor- 
mance is based on frequency instability as the basic phenomenon. CCIR 
Recommendation 538 therefore states the well known polynomial form of 
the spectral density of random frequency departures. It is now 
possible to estimate a probable time Interval error if the model para- 
meters of the clock are known. This error depends on the initial 
synchronization error, the frequency drift, the initial frequency 
setting error and the random frequency Instability. 

With established relations between the frequency instability of clocks 
and the timing properties required for the system, the relevant param- 
eters determining the timing properties can be Identified. 

In a network these parameters are a) the characteristics of the clocks 
located at the nodes and b) the spectral densities of the delay 
variations of the links. One Important optimization problem is the 
design of remote control systems for clocks from a master clock over 
links affected by delay variations and noise. The result depends on 
the bandwidth of the control loop and the free running stability of 
the controlled clock. 
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1. INTRODUCTION 


The purpose of this paper is to present some definitions, relations 
and practical data for the planning and design of network timing in 
digital communications systems. It is based on information which in 
most cases has already been published in various forms and also provi- 
des an update on the current work on recommended performance standards 
which is going on in the International Telephone and Telegraph Con- 
sultative Committee (CCITT) and the International Consultative Com- 
mltee on Radlocommunlcatlons (CCIR) of the International 
Telecommunications Union (ITU). 

There is a need for clarification of concepts and unification of 
language among the people working on new devices and systems related 
to the field of PTTI Applications. 

The system planner and designer needs well established theoretical 
concepts. In addition, he needs a good feeling for the orders of 
magnitude for the relevant parameters. The data presented in this 
paper are mainly illustrative and not intended to promote a particular 
design or product. 

In section 2 clock parameters are summarized. Section 3 deals with the 
concept of Time Interval Error and section 4 with transmission delay 
variations l.e. the degradation of timing information occurring on 
communication links. In section 5 the performance of a remote 
controlled clock is discussed by means of numerical examples. 

2. CLOCK PARAMETERS 

The performance of clocks is described by the following set of 
parameters. 

- Environmental Influences: - temperature 

- shock 

- acceleration 

- humidity 

- Frequency drift or aging 

- Random frequency instability 

Environmental Influences limit the performance of free running clocks. 
In the case of frequency controlled clocks, the control range will 
have to be large enough in order to compensate for these effects. 
Environmental specifications therefore constitute an Important set of 
design parameters which often determine the choice of a system con- 
figuration and the requirements on the devices used. The values given 
in the numerical examples are for a mild environment typical for clvi- 
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lian applications (Temperature +10 to +50“C, Humidity 30 to 80 percent 
relative, negligible shock and acceleration). Military or civilian 
mobile and aeronautical applications require other sets of parameters. 

Assuming thus that we have environment under control, we can con- 
centrate on the other two parameters, whereby these specifications 
Include some residual environmental effects (e.g. the flicker 
frequency Instability In most operational clocks). 

Table 1 contains stability data on the following types of clocks 

cesium 

rubidium 

crystal, oven controlled, high stability types 

The aging coefficients are upper limits derived either from 
manufacturer's specifications or from operational experience. Four 
different types of crystal oscillators have been considered. 

No. 1 Is a 5 MHz commercial type In wide use today 
No. 2 Is a less expensive small oven controlled VCXO 
No. 3 Is a low noise 10 MHz VCXO 
No. 4 Is a projection of new designs 

The random Instability coefficients are those of the well known 
theoretical model [1-6] recommended by the CCIR [7] and based on the 
spectral density of normalized random frequency departures: 

S (f) = h + h f"^ + h-f° + h-f^ + h-f^ (1) 

y -2 -1 0 12 

V-Vq 

with y(t) = — — 

0 

V being the actual "Instantaneous” value of the frequency, Vq the 
nominal value and f the ”Fourler"-frequency. The values In the table 
are computed from the time-domain data given by the manufacturers as 
the sqiiare root Oy(T) of the two-sample Allan-varlance , using the 
relation: 


a^(T) = h + h 2Zn2 + h, 

y -26 -1 0 2 t 


1 

+ h^ “Y7 <4.5 + 3£n (2TTf^T) - in2) h^ — ^ 


4tt t 


4tt t 


( 2 ) 


where f^ la the upper cutoff frequency of the measurement system used 
to measure a^(T), and T the sampling time. 
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The figures given are all conservative upper limits, usually 
guaranteed by the manufacturers and not so-called "typical” values 
which often tend to be optimistic. Only the crystal oscillator example 
No. 4 constitutes a projection not yet confirmed by large scale opera- 
tional experience, but prototype results are available [8]. 

The values missing from the table have not been Included because they 
are either not significant for the applications discussed In this 
paper (short term, high f fluctuations for Cesium and Rubidium) or 
small compared to the other terms (hQ for crystal oscillators). 

3. TIME INTERVAL ERROR 


Timing Is an essential function In any digital communications system. 
Therefore, International standards defining the minimum performance of 
the clocks used for network timing are required In order to achieve 
world-wide digital communication. The CCITT Study Group XVIII working 
on standardization of digital networks has Issued a Recommendation 
G811 In 1976 (9). A revised edition will be published In 1981 after 
the recent 6th General Assembly (Geneva 10-24 November 1980). This 
document specifies the minimum performance of clocks suitable for ple- 
slochronous operation of International, l.e. border-crossing digital 
links. Pleslochronous operation between networks means that each net- 
work Is controlled by a master or reference clock which Is free- 
running with respect to that of the neighbouring network. This 
requires close tolerance frequency adjustment and high stability In 
order to keep the rate of occurrence of slips below one In 70 days In 
any 64 kblt/s channel. 

The long term normalized frequency departure must therefore be less 
than 1 part In 10^^ In normal operation. The probability of degrada- 
tion to 1 part In 10^ shall be less than 10~5 and that of unavailabi- 
lity less than 10“^. 

The clock performance limits are defined as a maximum allowable Time 
Interval Error (TIE): 

TIE(t) = AT(t+T) - AT(t) (3) 

where AT Is the time difference between the actual timing signal and 
an Ideal timing signal. Using a more familiar notation [2, 3, 4, 6, 7] 
we have for any time Interval T: 


TIE(t ) = x(t +T) - x(t ) = / y(t)dt (4) 

^0 
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For a clock which has been adjusted to the nominal frequency and 
synchronized at t“0, the probable time interval error Cfjj(t) can be 
estimated. is the standard deviation of the time departures of 

an ensemble of identical clocks having random instabilities described 
by the model of equation 1. The following relation has been 
demonstrated by means of computer simulations [10]; 

= t (a^ + Oy (T.t))^/2 (5) 

OyQ is the standard deviation of the initial frequency adjustment and 
(TyVT) the two sample standard deviation (square root of equation 2) 
describing the random frequency instability. It is assumed that the 
parameters characterizing the clock do not change with time and that 
the initial adjustment error and the subsequent random frequency fluc- 
tuations are statistically Independent. In addition to the random 
instability and the initial adjustment error, frequency drift has to 
be considered and this leads to a combined systematic and random 
estimate of the TIE 

TIE = ^ + t (a^ + (T=t))^^^ (6) 

est. 2 '■ Yq y '' 

Eventually, significant cyclic environmental effects will have to be 
included, too. The method of estimating the TIE using equation 6 has 
been discussed in the CCIR Study Group 7 and Included in a draft 
Report [11] . 

The allowed limit of the TIE recommended in CCITT G811 is shown in 
Fig. 1. 

On the short term we see different curves which depend on the bit rate 
of the communications channel considered. This is due to the fact that 
the limit defined in this range is 100 t + 1/8 unit interval, 
expressed in nanoseconds where the unit Interval is the inverse of th< 
bit rate. 

In a real situation the random component of the time departure deter- 
mines a finite probability of violation which can be computed for the 
model described above. A recommended value will have to be defined 
after further study. 

4. TRANSMISSION DELAY VARIATIONS 


The systematic and random departures of the timing Instants from the 
ideal periodic timing waveform are designated by the terms of "jitter" 
and "wander". In the lltterature we find these departures referenced 
either to the signal period (phase jitter) or to a reference time 
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scale (timing jitter). The latter form is to be preferred since it 
does not depend on the particular bit rate. 

The distinction between jitter and wander is at present not very pre- 
cisely defined. One could treat it as a matter of Fourier Frequency 
since most specifications of jitter spectral density and jitter 
transfer functions stop at the low end around about 20 Hz. 

In a recent update on synchronization [12] one finds some examples 
extending the frequency domain to quite lower values for spectral den- 
sities having the character of white phase (or timing) noise. 

The word "wander" appears in some CCITT documents as a designation for 
slow systematic and random phase or time fluctuations. 

The time fluctuations observed at the receiving end of a transmission 
link are the sum of the fluctuations originating in the clock at the 
transmitting end and the path delay fluctuations [13]. The planning 
and design of a network timing system requires quantitative knowledge 
of both phenomena. As summarized in sections 2 and 3 we have a fairly 
complete picture on the modeling of clock performance. Reliable data 
on path delay variations are more difficult to obtain. This may be due 
to the large variety of transmission systems in use or under develop- 
ment and to the fact that long term path delay measurements on opera- 
tional transmission systems are much more difficult to organize than 
laboratory measurements on oscillators. Therefore, the data given 
below are neither complete nor definitive but rather an illustration 
of orders of magnitude to be expected and subject to some revision. 

The best known example of jitter accumulation on long transmission 
links is that occurring on the lowest hierarchy PCM links with a 
bltrate of 1544 kBit/s (USA and Japan) or 2048 kblt/s (Europe). The 
basic theory has been established many years ago [14] and recent 
measurements known to the author [15, 16, 18] appear to confirm the 
theoretical predictions. The spectral density of the timing jitter has 
the general approximate form shown in Fig. 2. This figure represents 
an envelope drawn above the spectrum which is more complex and has 
maxima and minima depending on the digital word structure or pattern 
of the signal. Using the envelope for design purposes is a conser- 
vative approach. The timing jitter spectral density envelope can be 
described by 


S^if) - l»20o(^ fg) 

where h2oo arid functions of the number N of cascaded regenera- 

tors and the bandwidth or Q of the regenerator timing circuit. 
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Fig. 3 shows this dependence on N for two examples: passive LC circuit 
regenerators with Q * 80 and PLL-type regenerators with Q “ 500. 

We see that for Fourier frequencies below fg, the PCM-llne jitter has 
the character of white timing (or phase) noise. 

For random, zero mean normally distributed test signals this property 
appears to be true for arbitrarily low Fourier frequencies. For real 
voice and data traffic the statistics can be different. However, 
modern coding techniques tend to approach random or pseudo-random 
signal statistics. Similar descriptions of transmission delay sta- 
tistics can be developed for all kinds of digital transmission 
systems. A complete review Is beyond the scope of this paper. 

Systematic slow delay variations are mainly caused by environmental 
Influences on the transmission medium. In the case of satellite 
systems, the dynamics of the system configuration [orbit parameters] 
play a major role. 

Temperature coefficients of metallic and optical fibre cables are sum- 
marized In Table 2. These figures are orders of magnitude and may vary 
depending on the cable type and manufacturer. A published reference 
[17] has been found only for the optical fibers, the other data have 
been gathered from various unpublished reports. The temperature 
variations show a dally and yearly cycle which depends on the climatic 
conditions. In temperate regions such as Europe an average temperature 
of +10“C, a peak-to-peak range of 20®C (year) and 2“C (day) can be 
expected . 

For transmission over a geostationary satellite, a 24 hour cyclic 
variation of the transmission delay Is due to Imperfect orbital para- 
meters (eccentricity. Inclination). Taking Into account the necessity 
of Improved stationkeeping due to the Increasing population on the 
geostationary orbit, a peak-to-peak dally delay variation of 600 
microseconds should not be exceeded. Long term drift and orbit reposi- 
tioning will cause additional variations [20]. Terrestrial microwave 
links over line of sight paths are known to be very stable, probably 
less than 5 nanoseconds variation due to the radio path Itself. 

A mostly unknown parameter Is the delay variation due to environmental 
conditions acting on the equipments used In the various transmission 
systems. Much work remains to be done and It Is surprising to note 
that J.R. Pierce's plea for Information on this matter In 1969 [19, p. 
629] went almost unnoticed. Whereas theoretical work Is abundant and 
sometimes elaborate In the published literature [12]., reliable 
experimental data on transmission delay variations still remains 
scarce. It Is therefore difficult to see which and what part of the 
many theoretical approaches are really relevant for operating system 
design. 
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For a given transmission system, the quantitative description of delay 
variations must comprise both random and cyclical components. The ran- 
dom component Is conveniently represented In the form of a timing 
fluctuation spectral density Sjj(f). It Is at present not known if low 
frequency divergent terms (flicker or random walk) are significant at 
very low Fourier Frequencies. If such effects exist, they are probably 
small compared to the effects caused by temperature cycles. Accumu- 
lated statistics on extended systems subject to many temperature 
cycles with different amplitudes and phases could result in a fllcker- 
llke spectrum having significant density down to frequencies of one 
cycle per year. There Is however a fundamental difference between 
transmission delay and clock time variations: The electrical length of 
the transmission path Is always finite and Its variations remain 
within physically limited bounds, whereas the time departure of any 
real clock grows without bounds. 

5. REMOTE CONTROLLED CLOCK 


The problem of remote control of the timing properties of a clock over 
a transmission link arises whatever the type of system organization 
chosen In practice, l.e. mutual synchronization, hierarchical master- 
slave (HMS) or any of possible combinations, refinements and self- 
organization schemes [21, 22]. The basic system elements are always a 
reference clock, a transmission link and another clock the time and 
frequency of which Is to be controlled by some means In order to 
obtain the fulfillment of a specification such as CCITT Rec. G811. The 
reference might be Itself a controlled clock or a free running master. 
The examples discussed In this section are on a quite elementary level 
compared to those discussed In recent publications [23-26], but 
Include spectral density models of the transmission time jitter and of 
oscillator Instability based on real data. 

In most systems currently In operation and under development, the 
hierarchical master-slave (HMS) type of operation Is applied. HMS is 
straightforward and fits well Into the structure of existing public 
networks at the local and regional level. Compliance with CCITT Recom- 
mendation G811 requires that master clocks are accurate In frequency 
with an uncertainty of less than + 1 part In 10^ and this implies 
reference to the frequency of UTC. Cesium clocks operated In a reaso- 
nable environment are accurate and stable enough to satisfy the 
requirements with only a minimum of surveillance and very few, If any 
readjustments. Rubidium and crystal clocks require Initial frequency 
setting and subsequent frequency control from an external reference to 
compensate for the inherent frequency drift (aging). 

Figures 4 and 5 show the time-domain frequency instability Oy(T) of 
typical cesium, rubidium and crystal clocks based on the data of Table 

1 . 
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Fig. 6 shows the estimated time Interval error for some clocks (see 
table 1) operated In the free running mode. It Is thereby assumed that 
each of these clocks has been adjusted In frequency and synchronized 
at t"0. 

The uncertainty OyQ of the Initial frequency adjustment for each 
type of clock Is assumed to be minimized, l.e. the averaging time for 
the measurement used for the readjustment Is chosen so that the uncer^ 
talnty corresponds to the minimum, flat portion ("flicker-floor") of 
the ay(T) diagram of Fig. 4 and 5 respectively. In some cases, noise 
In the process of measurement through a transmission link may lead to 
a higher value of OyQ and this must be taken Into account In such 
situations. Figure 6 gives an Idea of the possibilities of manual 
readjustment based on measurements. It allows to determine how long a 
lower hierarchy clock can be operated In the free-running mode within 
the specified limits of G811 after Interruption of the link to the 
master clock. It Is assumed that the frequency control system memori- 
zes the last valid setting of the oscillator control before the 
Interruption [27, 28]. It turns out that a cesium clock can be 
operated Indefinitely, a rubidium clock (a * 1*10“H per month) for 52 
days and crystal oscillators No. 1 (a ■ 3*10~H per day) for 1.6 days 
and No. 4 (a ” 5*10“12 pgj- day) for 5.2 days. 

To Illustrate further the operation of a remote controlled clock. 
Figure 7 shows a block diagram comprising a cesium master clock, a PCM 
link as mentioned In section 4 and a crystal oscillator controlled by 
the most elementary form of a phase-locked-loop, l.e. a first order 
PLL having a single pole RC loop filter. The jitter transfer analysis 
of the PLL Is made using the phase-time formalism described In [4, p. 
188 ff]. 

Table 3 shows the formulas used In the computations of Oy^(T) 
describing the random frequency Instability of the slave oscillator. 
Then the estimated time Interval error of the slave oscillator Is com- 
puted using the relation: 


TIE. (t) 
4 


t 0 

^4 


(T=t) 



( 8 ) 


The first term Is similar to equation 5 and the second term Is due to 
the fact that we use a first order PLL, a Is the aging coefficient of 
the slave oscillator. If the slave oscillator Is to be controlled over 
a noisy link, a low value of ("loose control") Is desired but then 
the second term shows the possible limit of Improvement obtainable for 
a slave oscillator having a given aging coefficient. 

The Input data for the numerical examples are given In the Tables 4 
and 1. The figures which follow show a selection of some typical 
cases . 
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Figure 8 shows the spectral density Sy^(f) of the controlled oscilla- 
tor output for the 4 slave oscillators considered, both values of 
shown in Table 4 and with and without the PCM line jitter. 

Figure 9 shows the corresponding time domain stability ay^(x). 

Figure 10 shows the TIE estimate obtained by means of equation 8. 

It is to be emphasized that these results do not constitute a design 
proposal. They are shown to Illustrate the relation between reference 
Instability, line timing jitter, control loop parameters and slave 
oscillator Instability using a simple PLL-Control System. 

In a practical design, step frequency control with a memory and pro- 
vided with a fast acquisition mode as described In [27, 28] or perhaps 
higher order or microprocessor controlled adaptive systems are or will 
be used. There Is not much difference between the loop transfer func- 
tion In our Illustrative example and those used In [27, 28], except 
for the larger computation effort required for the analysis of these 
more sophisticated designs. 

The difference between the computed TIE estimate and the design limit 
Imposed by Rec. G811 constitutes a margin In which the systematic 
cyclic delay variations remain to be Included. The allowed amplitude 
of these variations depends on the probability of violation allowed In 
the design. 

In the case of directed control HMS, the type of transmission medium 
limits the allowable distance between master and slave clocks. Double- 
ended timing [22] can In principle reduce the effect of slow delay 
variations but the delays Involved In the back and forth transmission 
of timing Information may change the transient response of the control 
system and careful analysis Is required In order to ascertain the sta- 
bility of the control. 

6. CONCLUSIONS 


The basic parameters relevant to the design of network timing systems 
describe the random and systematic time departures of the system 
elements, l.e. master (or reference) clocks, transmission links and 
other clocks controlled over the links. Using the definitions and 
notations recommended In CCIR and CCITT texts, the quantitative rela- 
tions between these parameters have been established and Illustrated 
by means of numerical examples based on available measured data. The 
examples have been limited to a simple PLL-control system but the ana- 
lysis can eventually be applied to more sophisticated systems at the 
cost of Increased computational effort. 

The author gratefully acknowledges the contribution of Kurt Hllty %dio 
did the programming for the computations of the numerical examples . 
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TABLE 1 


Type of clock Aging Random instability coefficients 




M(8) 

ho 

hi(s“l) 

h2(8-2) 

Cesium 

<5*10~l^/year 

2,885‘10”26 

<M 

1 

0 

• 

00 

- 

- 

Rubidium 

<1*10“^^ /month 

1,803-10"25 

5*10-23 

- 

- 

Crystal 1 

< 3 . 10-11 /day 

7,2*10-25 

- 

9,4*10-25 

3,13*10-26 

2 

<5* 10”^®/day 

2,89*10-22 

- 

- 

8,42*10-25 

3 

< 5 • 10“^®/day 

7,23*10-23 

- 

7,88*10-25 

7,88*10-27 

4 

< 5*10“12/day 

2,89*10-26 

_ 

9,4*10-26 

3,13*10-27 


TABLE 2 


Type of cable 

Temperature coefficient 
of delay In ns /km x °C 

Symmetrical pair, wire diameter 
0.8 mm paper Insulation 

3. ..5 

Symmetrical pair, wire diameter 
0.8 mm polyethylene Insulation 

0.3. ..0.75 

Coaxial cable polyethylene/alr 
Insulation 1,2/4, 4 mm 2, 6/9, 5 mm 

0.03. ..0.075 

Optical fibre S102 core [17] 

0.035 
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TABLE 3 


Input data; 
Source : 



© © © 



Line 


H (jw) 
X 




Derivations: S 


S 

Line 



Output: 



S + S 
^1 ^Line 



+ S 


40 



Transformations : 



2 2 
4tt f 




ITT 




. 4 
sin u 

2 

u 


du 


with u = TTfT 


General Rela- 
tions used: 


S = 4TT^f^ S 
Yi 
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TABLE 3 (CONTINUED) 


PLL-Data : 


DC Loop Gain K [s ^1 

X 

Single pole RC Filter 
Time Constant T^ [s] 


Damping ratio: ^ = — — — 

2/k T, 

X 1 


natural frequency f = — 
^ n 27T 


Loop transfer 
function: 



1 

2 2 2 2 
(i-f2 ) + 4Q r 


ll - hJ2 = |h^|2 


n = 


with 


_f 

f 


TABLE 4 


—26 —1 —21 
S = 2.885x10 f + 8.0x10 

^1 

-16 N = 552 

S = 2.49 • 10 • max. Jitter . 

X, . Q * oO 

Line 


PLL parameters 


1) K = 10 ^ C = — f = 2.3x10 Hz 

2) K = 0.1 s“^ ^ f = 2.3xl0~^ Hz 

"" ^/2 ^ 


(Crystal Oscillator Parameters see Table 1) 


Measurement system cutoff frequency 
for time-domain instability a (T) 

y 

f = 1000 Hz for all cases. 

H 
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Time Interval Error Limits, CCITT Rec. G811 (Draft Revision 1980) 












Fig. 4 


log «y(T) 


Cesium and Rubidium clock instability 


Fig. 5 



Crystal clock instability 

log Sy (T) 
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Fig. 6 TIE for free running clocks 
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Fig. 7 


REMOTE CONTROLLED CLOCK 
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Fig. 8 b) Controlled oscillator frequency domain instability, slave oscillator No. 2 
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Fig. 9 d) 


Controlled oscillator time-domain instability, slave oscillator No. 4 
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QUESTIONS AND ANSWERS 


MR. GEORGE PRICE: 

I have drawn the conclusion that one ought to spend a lot of money 
to get a very, very good high stability crystal oscillator for a 
digital communications system, because the requirements of CCIR or 
CCITT are such that we need that. Is that a wrong conclusion? 

DR. KARTASCHOFF: 

It is a right conclusion. You can use crystal, or rubidium, or 
cesium, you can also use very clean links and a lot of crystals. 

There are many, many possible solutions. Of course, we would 
welcome having a little bit better crystals. 

MR. PRICE: 

I could get by with a very cheap oscillator or I could get by with 
a real expensive one, and whether I spent the money to get a very 
expensive one would probably depend upon what? — accuracy require- 
ments, or slip requirements for the communications? I just don't 
know why I would spend the money to get a more expensive oscillator, 
is what I am saying. Is there a communication efficiency advantage 
in going to a very, very highly stabilized BVA type resonator that 
you recommended here a bit ago? 

DR. KARTASCHOFF: 

Well, I do not recommend it, I just have seen that it is a possi- 
bility that probably will come on the market, and of course it 
makes things easier, but the communications system will also work 
with the present oscillators. 
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